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ABSTRACT 


Theory for aeolian tone intensity created by rigid 
circular cylinders has been found to agree quite well with 
experimental results, his investigation is to exerimentally 
explore the wake of a cylinder shedding vertices at a frequency 


that is equal 
the cylinder. 
areas; 


1. 


2. 


to one of the natural vibrational frecuencies of 


Investigation was carried out in the following 


Patewal spacial correlation of vortex 
shedding, 


distance behind cylinder and valuc of 
RUS velocity at assumed position wnere 
WOmbIGesesmart. (GO) Told aii a 


observation of the phase between velocity 
Signals at two points along a cylinder. 
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INTRODUCTION 


The study of fluid flow past circular cylinders has been 
in progress since the latter part of the 19th century when Strouhal 
first showed that the vibrations causing the Sound were transverse 
to the direction of the wind, The early interest in tne phenomena 
was mainly concerned with the vortex shedding frequency as a 
function of Reynolds number, 

Until recently, very little work has been done concerning 
the intensity of vortex shedding sounds known as aeolian tones, 


(2) 


Pitts and Etkin, Korbacher and Keefe have derived 


theoretical expressions for the radiated sound intensity by starting 


(3) 


with Lighthills basic theory of aerodynamic sound as had already 


been extended to rigid bodies by ae 

Phillips used the product of 2 - to account for the 
fact that the force may not be coherent over the entire length of 
the cylinder. The use of the product 2& h is valid if the body is 


long compared to the coherent length, His resulting equation for 


sound intensity is: 


Po ae ea 5 
ab ea een oie S 2 Cc 2 P U 5 Q Q S 5 n =o 
16r | 


(1) 





), 


indicated time averages 


p = free stream density 
r,9 = polar coordinates (9 being the upstream direction). 
Nd 
5 i ne Strouhal number 
O 


d = cylinder diameter 
N = vortex shedding frequency 
U = free stream velocity 


Cy meeliic Coetricrent 


M = Mach number 

& = Jength of cylinder shedding vortices 
h = spanwise coherent length of vortices 
I = acouStic intensity 


Etkin, Korbacher, and Keefe derived a more general equation that 
considers any length of coherent force, 

In order to allow for longer lengths of coherent foree, 
the centroid of the area under the correlation coefficient nust 
be included, The effective length that is substituted for 27 is 
he (2-7) rather than & he 

Y = centroid of the area under the cross correlation 
coefficient curve, yy becomes small as the length of coherent 


Force becomes small, Therefore, equation (1) is apporached for 





small Le The more general equation from Etkin, et al, shows: 


I = ae 5? one Mm? ure (ny eedne 


(2) 
Experimental investigation of radiated sound intensity 
from circular cylinders when shedding vortex wakes has been made 


A@ (6) 


by Gerrar and Keefe » In general, these experimental re- 
sults confirm the theory for rigid bodies eSpecially in the 
lower Reynolds number region, 

The flow past cylinders that are nonrigid bodies has 
not been investigated sufficiently to determine what happens when 
the vortex shedding frequency and a natural frequency of the 
vortex producing body coincide, This thesis is a report of 
experimental investigation of certain properties of the wake 
behind a cylinder shedding vortices at a natural vibrational 
Frequency of the cylinder producing the vorticcs, 


The theoretical equation mentioned earlier for sound 


mrecnsity can be simplified for evaluation to: 


wre 4 


2 
Il & Cy Cas 
Cay 


which shows the three terws that may change for the situation 


of frequency coincidence Co, as and Y¥Yp Current work using the 





Same experimental set up is procceding to obtain experimental vajues 
for sound intensity, lift coefficient, and any "locking together" 

of the two frequencies, Bee expresses the opinion that the 
correlation length may increase preatly when the vortex shedding 
frequency and the vibration frequency coincide due to a coupling 

of fluid motions at different axial positions caused by the cylinder 
motion, This supposition was explored experimentally and the re 


sults are presented here, 
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CHAPTER IL 


PROCEDURE 


Ti,1l. Lateral Spacial Correlation 


RE FoF 


The aeolian tone test rig, which is pictured in 

Appendix A, was used to hold a cylinder of .154 inches diameter 
and 15 inches in length under very high tension to produce the 
desired natural frequency, It was desirable to have one of the 
lowey harmonic vibrational frequency coincide with the vortex 
shedding freauenctes in order to eliminate as many nodes jn the 
cylinder’s vibration as possible, The cylinder was put under as 
high tension as the test rig was physically capable of holding, 
The table below gives the natural frequencies of vibration and 
the air speed and Reynolds number required to cause vortices to 
shed at these frequencies, Only the wake of the 2nd, and 3rd, 


harmonic was explored, 


Ud 

. onal = p Et caren 

N (cps) U (fps) Re . 
2nd. harmonic 10390 62.5 ieee Wie oe 
3rd. harmonic OZ) 83 607 =x 10° 





Tne cross=correlation coefficient was computed in the 


following manner: 


Vives (an? - | — 4 AB 


2) nese ean 


me aN: re V Be 


R 








abs 


(4) 
Rap 64 = the cross<correlation coefficient as a function of 
spacial separation of the hot wire anemometers,. 

The cross-correlation coefficient is stated here to be 
only a function of the separation and not a function of where the 
separation is located along the cylinder. The investigation to 
feeermime at tite cociticient andjin turn, the correlation length 
is independent of location along the cylinder was not carried out 
and is a subject for future investigation, 

The values of R(z) were then plotted against z. See 
Fisures [4], [5], [6], and [7]. A mechanical integrator was used 
to determine the area under the R(€z) curve. This value is defined 
as one half the correlation length. 


oo 
iG, a 2 JS Ry, (4) dz 


om OO 


(5) 


A plot of ue as a function of Reynolds number is shown on Figure Fone 
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In this range of Reynolds numbers, the Strouhal number is practi- 
cally constant at (205, The cylinder consisted of a ,091" diameter 
steel wire over which was placed a Teflon tube that brought the 
overall diameter to .154 inches, 

This cylinder was placed in the open jet test section 
of the low turbulence wind tunnel in the N,1.2. Acoustics and 
Vibrations Laboratory. Figure [1A]shows the dimensions of the 
wind tunnel, It is more fully described in a report by fencer 

For the correlation measurements one of the hot wire 
sensing elements was placed in a stationary position while the 
other sensing element was held on a traverse mechanism that could 
be moved in a direction parallel to the cylinder, The two hot 
mares were oriented parallel to the cylinder and positioned at 
me edge of the wake where the signal is periodic, but not 
Eirbulents Tnhis position varied from speed to speed, but in 
general was about 1/2 inch back and 1/4 inch to one side of the 
center of the cylinder. See Figure [2]. The signals from these 
two Constant Temperature Anemometers (CTA) went to a correlation 
unit that produced a sum or difference oe Tne RMS output was 


recorded on a Graphic Level Recorder, Four outputs were recorded 
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To further check the correlation at coincidence of 
frequencies a dual beam oscilloscope was used to observe the 
phase relation of the signals from each hot wire, 

To ensure that the vortex shedding frequency is 
coincident with a natural frequency, a sine wave sipgnal from an 
oscillator, set on the desired natural frequency by an electronic 
counter, was displayed on the dual beam scope. The speed of the 
tunnel was then adjusted to cause the hot wire signal to match 
Eat {reduency of the oscillator, 

After setting the wind tunnel speed, the signals from 
both CTA's were displayed on the oscilloscope for observation of 


the phase between the two signals, 





fin?) OSittion of. liaximim RMS Velohity 
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To determine the poSition where it is assumed that the 


t ° ° 
one hot wire anemometer was used with a 


Mortrces start to toll up, 
two directional traverse mechanism, A search technique was used in 
the following manner, Starting very close behind the cylinder, the 
wake was traversed noting the position at which the maximum RMS 
velocity occurred. The value of RMS voltage and the DC voltage 
output of the anemometer was recorded at this location along with 
the position, The hot wire sensing device was then backed away 
from the wire a little further, The same traverse was made 
recording the same values at this point, This wes repeated until 
the maximum value of RMS voltase was bracketed, Then, the exact 
position was located by searching within the bracketed space. A 
typical plot of how the maximum RMS voltage varies as a function 
of distance behind the wire is presented as Figure [9]. The maxi- 
mim value is presumed to be the location of where the vortices 
Seent rolling wp’. 

A traverse behind the cylinder was performed at three 
speeds to have an idea of the properties of the wake, The results 


of RMS velocity and mean velocity are plotted and presented as 


Figure [10], [11], and [12]. 


CHARTER ETE 


DISCUSSION OF RESULTS 


iia. Lateral Spacial Correlation 


As stated in the introduction, the equation for the intensity 
Greacolian tones is proportional to lift coefficient, correlation length, 
and centroid of the area under the cross-correlation coefficient curve 


as thus! 


(3) 


The correlation bength a is the effective length of a 
coherent vortex, As the Reynolds number increases, the two 
dimensional qualities of the vortex shedding start to disappear. 

In other words, there is a span-wise variation in the phase of the 
vortex shedding, The average coherent length of shed vortices is 
used in calculating the acolian tone intensity, This average length 
is obtained by a procedure uSing two-point correlation of velocity 
along a cylinder as described in Procedure, 

Correlation of velocity variations rather than force were 
made for two reasons, First, it was impractical to put pressure 
taps on the cylinder duc to the small diameter, Secondly, because 
the cylinder has an induced motion the pressure tap readings would 
consist of a component due to the motion as well as the pressure 


caused by flow fluctuations, 


10 





3) 
ei Beroeeti. ) made similar cross~correlation measurements on 


a rigid cylinder, His data ranged in Reynolds number from about 
an to 4.5 x oe He found that the correlation length increased 
as Reynolds number increased, His value of de increasing from 
3.2d at Reynolds number toa to 6d at Reynolds number 4,5 x ily": 

Figure [8], shows correlation length decreasing with 
increasing Reynolds number from about 8d at Reynolds number of 
oe x 10° to 5d at Reynolds number of ll x to 

No apparent change in the trend of an was noted at points 
of coincidence of the vortex shedding frequency and a natural vibra- 
tional frequency of the cylinder, It was tnought that with the 
cylinder vibrating primarily in one mode that the coherent length 
would increase possibly approaching “ (n = mode number). Comma ae 


in fact, measured the correlation coefficient at low Reynolds numoer 


and found it to be of the following shape: 


R(z) 


(| 


ie 





In the measurements obtained during this investigation, 
the value of the correlation coefficient approached zero ionotoni- 
cally at all velocities including those at which the vortex shedding 
frequencies coincided with a natural frequency of the cylinder, 

More cross-correlation measurements should be made with 
particular emphasis on having the vortex shedding frequency exactly 
match the natural frequency and also more measurements at frequen- 
cies near the natural frequency. 

When observing the signals from both hot wires simultan- 
eously on a dual beam bectlloscone. with the hot wire probes 
adjacent to one another the signals remained inphase as would be 
expected, As the probes were moved apart parallel to the cylinder, 
the two signals did not remain inphase nor show a constant phase 
shift. Instead the two signals would move in and out of phase with 
me apparent pattern, 

To check for phase reversal when crossing nodal points 
the hot wires were placed on each side of a nodal point at a point 


midway between the nodal points as thus for the 3rd. harmonic, 


CTA 
ae A 
a oy 
Se itl ae 
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Again the signal moved continuously in and out of phase 
as previously described, 

In all these observations, on the dual beam scope, there 
didn't appear to be any change of observed signal whether the 
vortex shedding frequency was in coincidence or not with a 


natural frequency. 
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III.2. Position of Maximum RMS Velocity 


The position of the maximum RMS velocity is assumed to be 
the point where the vortex is "rolling sows 

The maximum RMS voltage referred to is the unfiltered hot 
wire anemometer signal. The signal at the position of maximum RMS 
velocity contained random frequencies together with the vortex 
shedding frequency, A frequency spectrum was taken of the hot 
wire anemometer signal, The background noise was at least 15db 
lower at all times than the signal at the vortex shedding frequency, 
Therefore, the addition to the broad band sipnal was less than 
feecan Or .O0l volts,. 

As shown in Fipure [9] the magnitude of the fluctuating 
velocity builds up linearly as the hot wire is moved downstream 
to a maxinum then falls off showing the viscous diffusion of the 
wert Lei ty. 

Birkhort ‘?? assumes that the wake acts as an oscillating 
flat plate attached to the cylinder providing alternating lift as 
the wake swings from side to side, He develops a satisfactory 
Strouhal number using this assumption, If it is assumed that the 


1 


point of inaxinum RMS voltage is the point where the vortex starts 


"rolling up" and that Birkhoff£'s evaluation of the wake acting as 


emeoscillator providing alternating lift, then the distance 


downstream to the point of maximum RUS voltaece and also the 
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distance off to one side of centerline where the maximum RMS voltage 
occurs would affect the lift produced by the wake, 

The distance downstream to the point of maximum RIMS 
voltage is seen from Figure [13] to decrease linearly in the 
veynolds number range considered, This follows the same trend as 


reported in Gages! 


from data taken by Schiller and Linke, 

Again we see no deviation from this trend when the vortex 
shedding frequency and a nautral frequency coincide, 

Figures (10, [11], and [12] show the mean velocity and 
RMS velocity of the wake at three different speeds with one of the 
situations Her aeene vortex shedding frequency coincide with the 
3rd, harmonic of the natural frequency, The results show that 
with heynolds number increasing the maximum difference in mean 
velocity increases indicating an increase in vortex strength, 

The maximum RMS velocity occurs at the point where the 
mean Velocity is halfway between the maximum and minimum inean 
velocity. This is expected because this is the area of maximum 
shear flow. Tne slope of the mean velocity increases in the 
region of maximum shear flow with an increase in Reynolds NUMDEX « 

In cenclusion, there did not appear to be any abnormal- 
ties in the properties of the wake when coincidence of frequencies 


eccurred, 
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APPENDED 
EQUIPMENT USED FOR CORRELATION LENGTH MEASUREMENTS 
Thermo-Systems, Inc, 

Model LO10A Constant Temperature Anemometer 

Model OSB “Corre lator 

Model 1031-2 Monitor and Power Supply with 

(2) tedel 1035 Linearizer modules, 

The sensing device used was a Tungsten Hot Wire .00015" 
diameter and .05" long. 

The output of the Constant Temperature Ancmometer was 
led into the linearizer modules so that given fluctuations in 
velocity would give the same voltage fluctuations for any mean 
velocity. The linearized output was sent to the correlator unit 
that summed and differenced the two signals:also the separate 
signals were obtained from the output of this unit. The output 
from the correlator was sent to the Bruel and Kjaer Frequency 
Analyser connected to the Graphic Level Recorder on which the 
meemeals were recorded in db level, Sce Figure {3}, The output 
was continuously observed on an oscilloscope. This was converted 
to RMS voltage and used in the cross~correlation coefficient 
equation. The magnitude of fluctuatiny velocity was not computed 


as the conversion factor cancels out of the R(€z) equation, 
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The free stream velocity of the wind tunnel was obtained 
by placing the hot wire sensing device in a free stream portion 
of the open test section and recording the output voltape from the 
Jinearizer to use on the previously obtained calibration curve, 
From this velocity, together with the cylinder diameter and 
kinematic velocity, the Reynolds number was computed, 

The method used to check for coincidence of frequencies 
was to run a frequency analysis of the output of an accelerometer 
mounted on a support of the wire. When the two frequencies 
coincided, the natural frequency would be reenforced to its 
highest value, A more accurate method of matching frequencies 


was devised described on page (8), 


2 





on 


© 


APPENDIX C 


RECOMMENDATIONS 


Take correlation coefficient measurements at different locations 
along the cylinder to determine if there is spacial homogenicty,. 
Work with a system that will allow a greater amplitude of 
cylinder vibration to determine how amplitude of vibration 
affects the correlation length. The cerrelation lengths should 
increase due to a coupling of mechanical and hydrodynamic forces, 
This could be accomplished by increasing the diameter of the 
cylinder without increasing the mass very much, This also 
would allow a coincidence of frequency at the lst, harmonic, 
Take more correlation measurement Jn close proximity to the 
point of coincidence of frequencies to determine if there is 
any deviation £rom the L vs, Reynolds number trend observed 

I GtioLleres1onc Ol 1 reqmeney. 

Take correlation measurement data in the same Reynolds number 
range that el Baroudi did to see what changes in trends or 


values are observed for a non-rigid cylinder, 
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